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Abstract. Vanadium dioxide VO2 and its metal-insulator transition at T=340K 
^ , continues to receive considerable interest. The question whether the physics of the 

insulating low-temperature phase is dominated by the Mott or the Peierls scenario, 
i.e. by correlation or band effects, is still under debate. A recent cluster dynamical 
mean field theory calculation [1! suggests a combination of both effects, characterizing 
' the transition as of a correlation assisted Peierls type. In this paper we present a 

detailed analysis of the excitation spectrum of the insulating Ml phase of VO2, based 
on this calculation. We implement a scheme to analytically continue self-energies 
^ ■ from Matsubara to real frequencies, and study the physics of the strong interactions, 

O . as well as the corresponding changes with respect to the density functional theory 

(LDA) band structure. We find that in the Ml phase life-time effects are rather 
negligeable, indeed allowing for an effective band structure description. A frequency 
^ ^ independent but orbital dependent potential, constructed as an approximation to the 

' full cluster dynamical mean field self-energy, turns out to satisfactory reproduce the 

, fully interacting one-particle spectrum, acting as a scissors operator which pushes the 

aig bonding and bands apart and, thus, opening the gap. 
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One of the key issues of condensed matter physics, both from the theoretical and 



^ ' the experimental point of view, is the understanding of the electronic structure of 
materials. With the advent of Density Functional Theory (DFT) within the Local 
Density Approximation (LDA), or its gradient extensions, the ab initio calculation of 
electronic properties of astonishingly many classes of compounds became feasible in a 
reliable fashion |2] . The manageable complexity of the material ever increases and even 
the verge of chemical precision comes into view. 

From a conceptual point of view, however, this might appear rather miraculous, 
as there are severe, well known, principle limitations. Indeed DFT is a theory of 
the ground state only. The interpretation of the Kohn-Sham eigenvalues of the 
effective one-particle system as excitation energies of the real N-particle system, has no 
theoretical justification, and only the striking compatibility with experiment promotes 
this interpretation, which is often tacitly assumed. 
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On the other hand, regarding the original domain of the theory, i.e. ground state 
properties, implementations of DFT inevitably have to resort to approximations, the 
exchange-correlation potential being unknown. This leads to the fact that for some 
even apparently simple materials, the calculated ground state is found to contradict 
experimental findings. This is the case when electronic correlations, i.e. genuine 
many-body effects, invalidate the band-picture which originates from an effective one- 
particle theory. The most prominent example is the Mott transition[3l H], where the 
correlations prevail upon the itineracy of the electrons, forcing them to localize. Among 
the materials where this mechanism is present are for example Cr-doped V2O3 or the 
parent compounds of the high Tc superconductors. A genuine signature of strong 
correlations, highly affecting the distribution of spectral weight, is the sensitivity to 
changes of external parameters, such as temperature or pressure, as is indeed observed 
in experiments such as transport, photoemission or optics. 

The formidable task is to unravel the behaviour of a system in-between these two 
limits : the effective free case (a problem that is diagonal in momentum space), and the 
atomic limit (that is diagonal in real space). In the latter the notion of bands having lost 
its sense, the quantity to discuss is the spectral function. Within a band-theory it will, 
for a given momentum, consist of Dirac distributions. With the increase of correlations 
these will broaden, due to scattering effects, inducing finite life-times, and loose weight 
that will reappear as incoherent features or satellites at different energies. 

In this work, we present for the specific example of the insulating phase of VO2, a 
study of how correlations modify the excitation energies of the system and how, if at 
all, this can be cast into an effective one-particle picture. 

To this effect, we implemented a scheme to analytically continue Matsubara self- 
energies to the real frequency axis, including for the first time the continuation of off- 
diagonal elements. This allows us to make a detailed analysis of recent cluster dynamical 
mean field calculations [1]. 

The paper is organized as follows : In section 2 we review the electronic structure 
of VO2 and recent cluster dynamical mean field calculations. In section 3 we apply 
the analytical continuation scheme (that we detail in Appendix A ) and discuss the 
real- frequency self-energies for the Ml phase of VO2. Section 4 presents a study of 
the correlation effects on the excitation spectrum and suggests an effective one-particle 
potential which corrects the LDA band structure. 



2. Electronic Structure of VO2 

2.1. Peierls versus Mott : A Reminder 

The phase diagram of vanadium dioxide VO2 has intrigued condensed matter physicists 
for decades[5l [6l[7l[8l[9l[l0l[IIl[12l[13l[14]. At high temperature VO2 is a (bad) metal 
and the crystal structure is of rutile type. At T=340K the compound undergoes a first- 
order metal-insulator transition (MIT)[15j and the structure becomes monoclinic (Ml 
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phase). The paramount characteristic of the low temperature structure is the pairing 
up of vanadium atoms, which form dimers along the c-axis that tilt out of the axis. 
Moreover the phase is found to exhibit no local moment behaviour. There is a long 
standing debate as to whether the MIT is primarily caused by the structural changes 
that double the unit cell (Peierls transition), or whether it is correlation effects that 
drive the system insulating (Mott-Hubbard transition). 

VO2 has a 3d^ configuration and hence the key focus lies on the d-orbitals. The 
cubic component of the crystal field resulting from the octahedral coordination of the 
vanadium atoms in the rutile phase splits the d-orbitals into an Cg and a t2g manifold, the 
latter being lower in energy. The total crystal field further splits the t2g into two Cg and 
a single aig orbital, which is directed along the c-axis. Both types of orbitals overlap, 
accounting for the metallic character of the rutile phase. In the Ml phase, where one 
has four formula units per unit-cell, the aig orbitals split due to the dimerization into 
bonding/antibonding states. Furthermore, because of an increased hybridization with 
the oxygen 2p orbitals owed to the tilting of the dimers, the are pushed up in energy 
relative to the former centre of the aig orbitals. Both these effects are working in favour 
of a gap-opening[6]. However, LDA calculations, though decreasing the density of states 
at the Fermi level as compared to the rutile phase, fail to capture the experimental gap 
(see figure Hj), which is only opened, once non-local correlation effects beyond the LDA 
are taken into account [U [161 flTj . 

2.2. Dynamical Mean Field Theory in a nutshell 

A theoretical model, that embodies the above mentioned struggle between localized and 
itinerant behaviour, or kinetic and potential energy, is the well-known Hubbard model. 
In the multi-band case its Hamiltonian reads 

H= 5^ <(k)cLck..+ E U^/^^i'^^^i''^' (1) 

k,«/',(T R,ll',crcr' 

where Hq' (k) is the one-particle dispersion, 1 being an orbital index, and U an on-site 
Coulomb repulsion that is here limited to density- density terms only, cj^^^, c^i^ are the 
usual fermionic creation and annihilation operators and n-m^ the corresponding number 
operators in real space. In a realistic context, the one-particle dispersion will be taken 
from e.g. an LDA calculation. Though the values of the interactions are often guided 
by constrained LDA|18j. RPA|19j or experimental estimates, they somewhat remain 
parameters of the technique. 

To solve this many-body problem in cases where the dominant correlation effects are 
local, a method of choice is the Dynamical Mean Field Theory (DMFT) and its realistic 
extension (LDA-I-DMFT) [20] . The tremendous merit of DMFT is that it is a non- 
perturbative technique that is able to accurately describe the local physics, independent 
of the interaction strength. The dynamical aspect allows for an accurate description of 
spectral weight transfer to incoherent features. 
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The basic idea of DMFT is to replace a lattice problem (or as in the realistic 
case, a problem defined within a localized basis set) by an effective one-site system, 
coupled to a bath and subject to a self-consistency condition, analogue to conventional 
Weiss mean field theory in statistical mechanics. Contrary to the latter however, the 
intervening mean field is energy-dependent, hence the notion of a dynamical MFT. 
For solving the DMFT equations, a variety of techniques are available. In the case of 
realistic multi-band calculations, a widely used technique is the Quantum Monte Carlo 
(QMC). A common feature of QMC implementations is that they work within the finite 
temperature Matsubara formalism, i.e. in imaginary time or frequency. Hence, an 
analytic continuation to the real frequency axis is needed. In Appendix A we describe 
such a procedure for the DMFT self-energy in the most general case, that for the first 
time includes the continuation of off-diagonal elements. The scheme thus gives access 
to the full orbital physics, as well as to non-local (Cluster DMFT) effects. 



2.3. Recent DMFT Results for VO2 

Recently, both the rutile and the Ml phase of VO2 were investigated within DMFTp!l[T71 
[T6] . Not surprisingly the insulating character of the monoclinic phase proved to be not 
capturable in the standard (single site) LDA-I-DMFT formalism without the unphysical 
creation of local moments. Instead, a cluster-extension of DMFT (CDMFT), where the 
solid is replaced not by a single impurity, but by a two-site cluster of vanadium atoms 
(surrounded by their oxygen octahedra) was needed[T]. It is of course very appealing 
from the physical point of view, to choose the vanadium dimers as the fundamental 
unit of the calculation. Herewith non-local inter-vanadium correlations are treated 
in a non-perturbative manner. The calculation in Ref [1] used an A^MTO-Wannier 
Hamiltonian[2T] where all orbitals other than the were downfolded, an on-site 
Coulomb repulsion U = 4.0eV, a Hund's coupling J = 0.68eV and T = 770K. For solving 
the DMFT equations, the Hirsch-Fye Quantum Monte Carlo (QMC) algorithm j22] was 
used. The insulating nature of the Ml phase is correctly described and the resulting 
local spectral functions of both phases compare favourably with photoemission and x-ray 
experiments [231 [24l[25l[26l[27l[28l[29l[30l[31]. The authors concluded that the non-local 
correlations effectively assist the Peierls-like transition, the dimerization leading to the 
formation of molecular singlets in the aig channel embedded in a bath. 



3. The DMFT Self-Energy on the Real Frequency Axis 
3.1. Orbital Structure of the Self- Energy 

In the present work, we have implemented an analytic continuation procedure for 



imaginary frequency DMFT self-energies (see Appendix A ), which enabled us to develop 
a more detailed analysis of the above calculation. With the real-frequency self-energy at 
hand, interesting quantities become accessible, among which figure momentum-resolved 
spectral functions and linear response functions, such as the optical conductivity (when 
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neglecting vertex-corrections). In the following we shall analyse the correlation induced 
modifications in the excitation spectrum of Ml VO2 and discuss whether the result can 
still be cast into an effective band-picture. 

With four formula units in the unit cell, the Ml VO2 self-energy matrix acquires 
the following form : 
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where the blocks Sjj are 3x3 matrices each, and S 



z = 1,2) are the local self-energy 
matrices for a given vanadium atom, and S12 and S21 are the inter- vanadium, intra- 
dimer blocks. Still, when given this additional freedom, the system is found to have in 
total only four non-equivalent self-energy elements. All other elements are zero within 
numerical precision. The self-energy for one dimer, i.e. the upper-left or down-right 2x2 
block in ([2]) comes out to be 
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i.e. the intra-dimer coupling is entirely due to the aig channel, as can be expected from 
the discussed symmetry of the system. 



3.2. Interpreting the Self-Energies : Correlation Effects in VO2 



Using the procedure detailed in [Appendix A[ we compute the DMFT self-energy on 



the real frequency axis. Figure [T] and [2] display our result. We have chosen to change 
the basis of the aig components and to show their bonding/antibonding combinations 
Tib/ab = ^aig ± ^aig-aig, in which the self-energy is diagonal in orbital spac^. 

Several important observations can be made. The suppression of spectral weight at 
the Fermi energy can be achieved by two distinct behaviours of the imaginary parts of 
the self-energy within the charge gap: Either they diverge (as is the case of the insulating 
phase of the one-band Hubbard model), or they have to vanish within the gap. In the 
present case of Ml VO2 we see, figure El that the self-energies are all regular and their 
imaginary parts vanish within the gap. Hence it is not a divergence of the effective mass 
that is responsible for the MIT. 

t We stress that this basis transformation can only be done after the DMFT calculation, as in this 
molecular picture the charges will no longer be centered on the atomic sites, and as a consequence the 
Hubbard-Hund Coulomb interaction terms would acquire an untractable form. 
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Figure 1. Real parts of the self- 
energy oig bonding/antibonding 
basis. The two stripes represent 
the position and bandwidth of the 
former LDA aig bonding (lu < 0), 
antibonding (w > 0) band. See 
14. 1.11 for a discussion. 
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Figure 2. Imaginary parts of 
the CDMFT self-energy of the 
Ml phase of VO2 in the aig 
bonding/antibonding basis. Note 
that all self-energy elements are 
regular at zero frequency. See 13.21 
for a discussion. 



Concerning the elements, both the real and the imaginary part, depend only 
weakly on frequency. Moreover the imaginary part is globally low in magnitude, 
signaling only minor life-time effects. This is a consequence of the low occupation 
of these orbitals (~16%). The interesting physics is carried by the aig components. 
They exhibit a rather selective dependence on frequency: The [antijbonding component 
has a quite flat ([positive] negative) real and an unimportant imaginary part in the 
[unjoccupied part of the spectrum. The d^-configuration, with its nearly filled bonding 
bands, and empty anti-bonding bands, thus leads again - but this time in a non-trivial 
way - to only small life-time effects. Despite the considerable on-site interactions, the 
complex structure of the self-energy thus suggests the existence of well-defined one- 
particle excitations. As we shall see below, this effective band structure is however 
strongly rearranged as compared to the former LDA energies, due to the differences in 
the regions of constant real-parts. 

4. Effective Band Structures : At the Rescue of the Band-Picture 
4.I. The one-particle Poles 

When neglecting the anti-hermitian parts of the self-energy, which - as we have seen 
above - are small for the bonding/antibonding bands, the excitation energies of the 
system are given by the poles of the one-particle Green's function : 

detK + /i-^o(k) -3?S(cUk)) =0 (4) 

In contrast to the DFT equations, the above is not an eigenvalue problem, due to 
the frequency dependence of the self-energy. As a first step we shall discuss a graphical 
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construction for finding solutions of (jlj). 

4-1.1. One-particle Excitations: A Craphical Construction. In the one-band case the 
problem reduces to cuk + — Ck — 3^S(ct;k) = 0. Solutions are the intersections of 
the real-part of the self-energy (minus /i) with a frequency stripe of slope one and 
a width corresponding to the dispersion of ey_. In the general Hamiltonian case this 
construction is no longer exact, but might still give qualitative insights. In the present 
case we expect this procedure to give a rather good description, when working in the 
bonding/antibonding basis and limiting ourselves to a discussion of the aig orbitals only. 
Indeed, if the orbitals were empty, this contruction would be exact. 

In figure [1] we indicate the positions and bandwidths of the aig bands within 
LDA. As in the bonding/antibonding basis the orbital pole equations approximately 
decouple, we show the LDA [anti] bonding bandwidth only in the region [u; > 0] a; < 0, 
where relevant intersections can be found. Below the Fermi-level, bonding pole solutions 
appear between -0.5 and -l.OeV, due to intersections with the bonding part of the self- 
energy (red) with the stripe representing the former LDA aig bonding band {uj < 0). 
As a result, the bonding orbital is considerably pushed downwards in energy. The aig 
anti-bonding poles, i.e. the intersection of the (blue) anti-bonding self-energy and the 
stripe at > are pushed beyond 2eV. Thus, already at this point, we suspect the 
spectral weight above 2eV as seen in the local spectral-function of Ref . [1] to stem from 
the aig anti-bonding excitation rather than from an upper Hubbard band. We will come 
back to this point in the next section. 

4.1.2. One-particle excitations : Poles of the Green's Function. To solve the quasi- 
particle equation (jlj), one could in practice scan through frequency and look for sign- 
changes of the determinant (jlj). However, as along various symmetry-lines, degenerate 
solutions exist, we had to introduce a threshold and solutions where identified as yielding 
an absolute value of the determinant below this value. Figure j3] shows our results. 

On general grounds, there are, at a given k-vector, at least as many poles as 
there were LDA bands (including possible degeneracies). However, the frequency 
dependence of the self-energy may produce additional solutions of (jlj) : In the particle- 
hole symmetric one-band Hubbard model, for example, there is one solution in the 
non- interacting limit, which simply is the one-particle dispersion, whereas in the atomic 
limit there are two poles at ±f//2. In between there is a regime in which three poles 
may occur, representing the famous three-peak structure in the spectral function. In our 
example of the Ml phase of VO2, we note that the number of one-particle DMFT poles 
equals the number of LDA-bands. In other words, the bands get shifted and reshaped, 
but the importance of incoherent features is rather small. This finding is far from trivial. 
For example in the metallic rutile phase of VO2P2] or V203[33j, additional poles occur, 
representing incoherent Hubbard bands. 

The most important changes as compared to the LDA is the opening of a gap 
between the bonding aig and the Cg bands. Besides, the anti-bonding a^ band is 
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Figure 3. Comparison between tlie one-particle poles (red) according to ([4]) and 
the bands using an optimized orbital dependent but frequency-independent effective 
potential obtained as an approximation to the full CDMFT self-energies (green) for 
various symmetry lines. See text for a discussion. 



shifted well beyond the ones, as already seen in the graphical discussion above. In 
a molecular picture, the bonding/antibonding splitting of the aig-orbitals, due to the 
dimerization, is tremendously enhanced by the non-local intra-dimer correlations that 
are caused by the local interactions. These results confirm the picture by Koethe et al 
who identify a high frequency XAS intensity peak as the aig anti-bonding band, while 
a small low-energy feature in the partial aig spectrum, both seen in the experiment and 
the CDMFT calculation results from e^-Oi^ hybridization [31] IT]. 

The neglected life-time effects will in our case broaden the spectra and shift the 
positions a little due to their frequency dependence, as can be seen in the (k-integrated) 
spectral function[T]. However they will not, as is e.g. the case in V203[33], completely 
suppress spectral weight for some given pole solutions. Hence, in the current example, 
the effective band structure gives a quite accurate guidance to the physics of the system. 

At first sight, it may seem surprising that despite strong local correlations - the 
imaginary parts of the on-site selfenergies reach values of the order of 3eV - well-defined 
one-particle excitations survive. Mathematically spoken, this is the consequence of 
a cancellation of local and inter-site self-energies resulting in nearly negligeable life- 
time effects for the [ ant i-] bonding bands at [positive] negative energies, respectively. 
Physically, enhanced inter-site hoppings allow the electrons to avoid the strong on- 
site repulsion by delocalizing over the dimer, and the resulting intra-dimer fluctuations 
reduce the net effects of the interaction [8]. 
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4-2. Tailoring of a one-particle Potential beyond the LDA : Scissors at Work 

As we have seen, the frequency dependence of the bonding/antibonding self-energy in 
regions where they yield solutions of the pole equation is rather unimportant. Hence 
the question arises as to which extent one could derive the effective band structure 
from a constant one-particle, albeit orbital-dependent, potential A. To this end we will 
approximate our dynamical self-energy in a static way, by its values at the former LDA 
band-centres for the e^, and by the values at the pole energies for the aig orbitals. This 
means we choose Ag^^ = 3?Se^^ (0.5el^) — fJ. = 0.48el^, Ae^_^ = 3?Se^^(0.5ey) — /i = 
OMeV for the components,' and A^ = 3?Sb(-0.75eV") - /i = -0.32eV, Aab = 
3fJSaf,(2.5el^) — /i = l.2eV for the aig bonding/antibonding components respectively. 
As a consequence, there will be a downshift of the aig bonding band (A;, < 0), which, 
together with an upshift of the bands (Ae^^^^ > 0) is responsible for the opening of 
the gap. Thus, these considerations put numbers to Goodenough's picture[6j. However, 
the effect caused by the re-hybridisations with the oxygen bands, as described in the 
introduction, is supplemented by the following : The correlations cause a depopulation 
of the Cg bands, in favour of the a^ bonding one, as thereby the electrons are effectively 
avoiding the on-site Coulomb interaction. Moreover, the anti-bonding aig band is pushed 
up in energy way beyond all the e^, tremendously enhancing the bonding/antibonding 
splitting. 

After transforming back into the original basis, we solve the eigenvalue problem for 
Ho(k) + A. The results are shown in figure [3] along with the exact poles from above. 
The shiftings of the bands, owed to the additional potential, as described above, are, 
besides minor details, in excellent agreement with the CDMFT poles. 

The above one-particle potential acts like a scissors-operator and is somewhat 
reminiscent of results obtained within the GW approximation [Mj. Indeed, within a 
model GW approach, a gap was found to open in the Ml phase[35]. Furthermore, 
preliminary results of a fully ab initio GW calculation suggest that the GW 
approximation is as a matter of fact suited to describe the insulating character in the 
present case. 

4-3. Linearizing the Self- Energy 

Though disposing of the fully frequency dependent DMFT self-energy, we find it 
instructive to compare with an approximate low-frequency expansion, one would 
employ if one had used a many-body technique that gives only information about the 
renormalization of low lying coherent excitations, such as e.g. slave-boson theories. The 
renormalized band structure is then given by the solutions cj^ of 

det (cjk - Z [Ho{k) + ^ 0) - /i]) = (5) 

where = 1 — di^QT,{iuj)\i^^o, and [i^SS] 3?S is the [anti-] hermitian part of the 
self-energy on the imaginary axis. This approach, exact at the Fermi-level and normally 
identified with the Fermi-liquid regime, may in principle be applicable in our case to the 



Effective Band Structure of Correlated Materials - The Case of VO2 



10 



insulating phase of VO2 as the self-energies do not diverge at low frequency, as is the 
case in the Mott-insulating phase of the one-band Hubbard model. Albeit, the quantity 
Z should, in this context, not be identified with the usual quasi-particle weight. 




Figure 4. Comparison of the LDA-bands (red) with the bands resulting from the 
hnearization scheme of (O (blue) 

As can be seen in figure HI a charge gap is found to open, thus correctly describing 
the insulating behaviour; however when comparing to the one-particle poles, we remark 
huge discrepancies: The overall bandwidth is too small, the bands too fiat and 
the bonding/ ant ibonding splitting lacks the tremendous increase found before. This 
is natural, since foremost the involved slopes of the self-energy at low-frequency 
overestimate the renormalizations towards the Fermi-level, whereas at high enough 
energy the self-energies are essentially constant. 

4.4- Conclusions 

In conclusion, we have presented a detailed analysis of the excitation spectrum of the 
Ml phase of VO2, with special emphasis on an effective band structure description. 
Indeed, we found life-time effects to be rather negligeable and the correlation effects 
to only shift and reshape the Kohn-Sham energies with respect to the starting LDA 
calculation. Furthermore, by approximating the dynamical self-energy, we were able 
to construct a frequency- independent, though orbital-dependent, one-particle potential, 
that, when added to the LDA reproduced the DMFT excitation energies to a quite 
satisfactory extent. We remark that due to the orbital-dependence of this potential, it 
cannot be viewed as an improved density functional. 

Hence, despite the undeniable presence of strong local correlations, the system 
retains the coherence of its excitations by means of intra-dimer fluctuations, to an 
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extent that the physics of the compound is indeed dominated by the Goodenough- 
Peierls-picture. The role of correlations consists in (i) pushing the aig anti-bonding 
band beyond the top of the e^, consistent with the experimental findings, and, more 
importantly, in (ii) enhancing the aig bonding - splitting due to an effectively reduced 
Coulomb repulsion in the aig bonding band, hence favouring the depopulation of the Cg 
bands. The latter results in the opening of the gap. Thus, as a matter of consequence 
it IS the correlations that are responsible for the insulating state, albeit they cause it in 
a rather specific fashion. 
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Appendix A. Analytical Continuation of the DMFT-Self-Energy 

Whereas for the continuation of the Green's function from the Matsubara formalism 
to real frequency a numerical technique, called the Maximum Entropy algorithm 
(MaxEnt)|36j has become a standard approach, the analytic continuation of the self- 
energy has only recently been envisioned in the context of realistic (i.e. multi-band) 
calculations [371 ES [39l HQ] . 

In this work we have implemented a scheme that is sufficiently general to work in 
the multiorbital cluster DMFT case, and which, in particular, allows for off-diagonal 
self-energy elements. The following sketch outlines our continuation procedure : 

QMC : Gair) 

I MaxEnt 
Au'iuj) = -1/7tQGu,{uj) 

I Kramers-Kronig 

Gu'iuj) = + /X - Hoik) - S(cu + ^0+)]-^,. (A.l) 

k 

I Root-finding 

Starting from the local Green's function in imaginary time, we perform the 
aforementioned MaxEnt algorithm, yielding the spectral function, which is proportional 
to the imaginary part of the real-frequency Green's function. The real part of this 
quantity is related to the former by Kramers-Kronig transformation. The by far most 
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difficult step in this scheme is the last one. Whereas in the one-band Hubbard model 
on the Bethe-lattice an exact expression for the self-energy in terms of the local Green's 
function exists (S = u + fx — t^G — 1/G), in the general multi-band case a multi- 
dimensional root-finding procedure ^41 J has to be employed to solve (lA.ip . In the present 
case, of four formula units in the cell, the Hamiltonian is a 12x12 matrix and the self- 
energy has the structure shown in section 13.11 

In the degenerate case, the k-sum in flA.l|) can be replaced by an integral over 
the density of states and one ends up with the simpler task of inverting the Hilbert 
transform. 

In the most general case, the starting Green's function will furthermore have off- 
diagonal elements in orbital space. As the standard implementation of the MaxEnt- 
algorithm in this context has, to our knowledge, so far been used for diagonal elements 
only, some comments are in order. The equation that is to be inverted reads 



Gu'ir) 



dujAii>{uj)Kf3{T,uj) 



(A.2) 



with Ki3{t,uj) being the usual MaxEnt kernel [36] . 

In contrast to the diagonal elements of the spectral function, which are normalized 
to one and are positive at all frequencies, off-diagonal elements have vanishing zeroth 
moment and hence change sign as a function of u. 

For mending this difference, the idea is simply to add a function i{uj) of zeroth 
moment one to const ■ A{uj), and calculate the corresponding shift in G(r). However 
this function has to be chosen such, that the resulting fictitious spectral function is non- 
negative for all frequencies. To this end, we construct it in the following way. Given the 
orthogonal orbital transformation A = WAU with 



/ 1 



1 

V2 



II' 
-1 



\ 



(A.3) 



1/ 

the 11-element {l^l') of the spectral function in the new basis reads 

~ All + A'l' 

All = ^ ± Ai, 



(A.4) 



and is thus per construction a proper diagonal element of the spectral function in the 
rotated basis, i.e. in particular it is positive for all u. Therefore / = l/2{Ai'i' + An) is 
a suitable shift for ^Aw in the original basis. 

Hence, after having continued the original diagonal elements, one constructs as 
input the function 

(A.5) 



Giiir) = {U^G{t)U)^^ = / dujAn{u)Kf,{T,LU 
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from which, using (lA.4p the off-diagonal spectral function is then deduced. We note 
that in the present case the above rotation has the same geometrical interpretation as 
the transformation into a bonding/ ant ibonding basis, as used for the a^g self-energies 
of VO2, but the construction is more general and does not require special symmetry 
properties of the system. 
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